EURONU-WP6-10-24 



MiniBooNE and LSND data: non-standard neutrino interactions 
in a (3+1) scheme versus (3+2) oscillations 

Evgeny Akhmedov 1, 2 '|j and Thomas Schwetz 1 '^ 

1 Max- Planck- Institut fiir Kernphysik, 
^ ; PO Box 103980, 69029 Heidelberg, Germany 

bJQ! 2 National Research Center "Kurchatov Institute" 

5 \ 123182 Moscow, Russia 

m 

^ Abstract 

Qh The recently observed event excess in MiniBooNE anti-neutrino data is in agreement with the 

LSND evidence for electron anti-neutrino appearance. We propose an explanation of these data in 
terms of a (3+1) scheme with a sterile neutrino including non-standard neutrino interactions (NSI) 
at neutrino production and detection. The interference between oscillations and NSI provides a 
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source for CP violation which we use to reconcile different results from neutrino and anti-neutrino 
data. Our best fit results imply NSI at the level of a few percent relative to the standard weak 
interaction, in agreement with current bounds. We compare the quality of the NSI fit to the one 



obtained within the (3+1) and (3+2) pure oscillation frameworks. We also briefly comment on 



using NSI (in an effective two-flavour framework) to address a possible difference in neutrino and 
anti- neutrino results from the MINOS experiment. 
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I. INTRODUCTION 



Recently the MiniBooNE collaboration announced updated results of their search for 
Vn —* v e transitions (l|. In the full energy range from 200 MeV to 3 GeV they find an 
excess of 43.2 ± 22.5 events over expected background (the error includes statistical and 
systematical uncertainties). In the oscillation-sensitive region of 475 MeV to 1250 MeV 
the background-only- hypothesis has a probability of only 0.5%Jl|. This result is consistent 
with the evidence for v^ — \ v e transitions reported by LSND [2] . Any explanation of these 
hints for v^ — > v e transitions at the scale of E/L ~ 1 eV 2 has to satisfy strong constraints 
from various experiments. First, no evidence for transitions has been found in MiniBooNE 
neutrino data above 475 MeV J3j. This suggests that CP (or even CPT) violation has to 
be invoked to reconcile neutrino and anti-neutrino data. Second, severe constraints exist 
for P« 0, | and „„ P„ @fl disappearance at this seale, which have to be respected by any 
explanation of the z/ M — > v e excesses. 

The standard approach to the LSND problem is to introduce one or more sterile neutrinos 
at the eV scale 10h12|. Adding one sterile neutrino one obtains the so-called (3+1) mass 
scheme. In this framework there is no CP violation at short baselines, and disappearance 
experiments strongly disfavour an explanation of the appearance si gnal s, see for example 13, 



14 



If two neutrino mass states at the eV scale are present [15|, [16| ((3+2) scheme), the 
possibility of CP violation opens up (l7j, which allows to reconcile LSND and MiniBooNE 
neutrino data 18j. However, constraints from disappearance data still impose a chal leng e to 
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the fit, and the overall improvement with respect to the (3+1) case is not significant 

In the following we will update the (3+2) results of 18] with respect to the latest data from 

MiniBooNE anti-neutrinos. 

Apart from sterile neutrino oscillations, various more exotic explanations of the LSND 



signal have been proposed, among them, neutrino decay 20 
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CPT violation 
violation of Lorentz symmetry |26H28||. quantum decoherence |2 
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neutrinos 
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mass- varying 
, shortcuts of sterile neutrinos in extra dimensions |33| or sterile neutrino 
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oscillations with a non-standard energy dependence 

In this work we attempt an explanation of the global data by departing from a (3+1) 
neutrino scheme through the inclusion of non-standard interactions (NSI) of neutrinos on top 
of the Standard Model weak interactions. Such new interactions may be induced by generic 
new physics beyond the Standard Model and can be neutral current (NC) like (involving 
only neutrinos in the lepton sector) or charged current (CC) like (involving a charged lepton 
and a neutrino). Model- independent bounds on such new interactions are at the level of 

and 
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few xl0~ 2 compared to the standard 4- Fermi interaction strength set by Gp, see 
references therein. An observation of NSI at that level would be a remarkable sign of new 
physics. The realization of NSI in terms of effective operators has been discussed recently 



m 



36, 
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Since the experiments considered here typically have rather short baselines (below 1 km), 



matter effects are very small and NSI affecting the propagation of neutrinos through matter 
will have a negligible impact. Therefore, we focus on CC-like NSI in the production and 



detection processes, see for example |38N41|. As pointed out in [39|, such NSI provide a 
new source of CP violation. Here we use this effect to reconcile neutrino and anti-neutrino 
data by exploring the interference of the NSI with sterile neutrino oscillations, with Am 2 of 
order 1 eV 2 . This allows us to keep NSI at the percent level. Previous attempts to explain 



the LSND excess by NSI without sterile neutrinos can be found in 42M44J . while sterile 



neutrinos with an exotic matter effect have been considered in (45|, |46 1 . 

In the following section [III we derive an effective parameterisation of the transition and 
survival probabilities for short-baseline experiments, identifying the particular combinations 
of mixing and NSI parameters relevant for the experiments. We show that the probabilities 
for the global short-baseline data in a general (3+1) NSI model depend only on 8 indepen- 
dent parameters. In this case we can make use of the fact that in LSND and KARMEN 
the neutrino production mechanism is different from all other experiments. Therefore these 
experiments decouple to some extent from the rest of the global data. On the other hand, 
neglecting NSI involving the charged muon, only 5 independent parameters remain, which 
are only two independent parameters more than in the case of standard (3+1) oscillations, 
one modulus and one phase. In this case the fit is more constrained, while still providing 
significant improvement with respect to the pure oscillation case. As we will see, this con- 
strained version requires only one single non-zero NSI parameter. In section IIHI we present 
the results of global fits within these two versions of the (3+1) NSI model, the constrained 
5 parameter version in section Ull Al and the general one in lHIBI Furthermore, in section HVl 
we present an updated fit in the (3+2) oscillation scheme and we compare the quality of 
the (3+1) NSI fits to the one of (3+2) oscillations. A summary and discussion follow in sec- 
tion [V] In appendix [A] we relax some of the assumptions made on NSI parameters (showing 
that they do not change our results significantly), and in appendix [B] we comment briefly 
on using NSI to address a possible difference in neutrino and anti-neutrino results from the 



MINOS experiment J47| 



II. NSI IN THE (3+1) MASS SCHEME 

A. The formalism 

We assume that, in addition to the standard charged current (CC) weak interactions, 
there exist non-standard CC-like interactions, whose Lagrangian can be parameterised at 
low energies as 

£ NS i = -2y/2G F J2 £ a^f P L,R^f')( J a PLl^) + h.c. . (1) 

a,/3 

Here G F is the Fermi constant, / and /' correspond to either quarks or leptons differing 
by one unit of electric charge and l a corresponds to a charged lepton (l a = e,/i,r). Pl(r) 



denotes the projection operator on left-handed (right-handed) fields. The particular chirality 
structure assumed in eq. ([[]) allows for interference of standard and non-standard processes. 
To simplify the notation we leave the dependence of e on L/R implicit. 

The interactions in eq. ([T]) contribute to CC processes of neutrino emission and absorption. 
In what follows, we will replace the superscript //' by X (X = S, D), where X = S stands 
for a concrete neutrino production process in the source and X = D stands for the neutrino 
detection process. In the presence of NSI the neutrino state \v%) produced or detected in 
a CC process along with a charged lepton l a is no longer the flavour eigenstate u Q , but a 
linear combination of flavour eigenstates: 



\v X 



> = ^(i^>+E e 5»i^>)- ( 2 ) 



A similar relation holds for anti- neutrinos, with e^p replaced by e^S. Note that in our 
notation the first subscript at the parameter e*o always refers to the charged lepton and the 
second to the neutrino. C* is a normalisation constant, which satisfies 

|Cf| 2 (l + 2Re^ Q + 5>*| 2 ) = l. (3) 

p 

We also assume that, in addition to the standard three neutrino flavour states v e , v^ and 
u T , there exists a fourth light neutrino, which necessarily has to be a sterile neutrino u s . As 
usual, the flavour states u e , u^, v T and v s are linear combinations of four mass eigenstates: 1 

!"«> = £ *>><>• ( 4 ) 

i 

We assume that the 4th mass eigenstate is separated from the other three by a mass gap of 
order Am| x ~ 1 eV 2 ((3+1) scheme). Alternatively, we will also consider a variant with two 
sterile neutrinos but no NSI ((3+2) scheme). 

The state \v%) can also be expressed in terms of the mass eigenstates: 



1^ 



') = c * E (^ + E e ^ u i) \**) = c * E (^ + £ i) u U»i) ■ (5) 

The probability amplitude for a neutrino u^ born in a production process in the neutrino 

D 



source to be detected as u? at a distance L from the source in a detection process is 



AapiL) = (uf\u s a (L)) =J2 F °i F £* e ~ iEiL > ( 6 ) 



1 In the presence of NSI given in eq. ([]}, the usual flavour states can only be defined at the level of the 
Lagrangian, as the states coupled to the Standard Model W boson. Therefore, the mixing matrix U is 
the unitary matrix which diagonalises the neutrino mass matrix in the basis where the Standard Model 
CC interaction (without new physics contributions) is diagonal. 



where Ei is the energy of the ith mass eigenstate neutrino and 

^[ = tf« £(*<* + <)*£■ (7) 

p 

Next, we note that for baselines and neutrino energies of interest, one can neglect all 
mass squared differences except those involving the fourth mass eigenstate. Denoting 

A^^IiL, (8) 

2E ' v ; 

we can rewrite eq. ([S]) in the generic form 

Aafi(L) = a a(3 (e~ iA - 1) + (5 aP , (9) 

where 

««/) = F a4 F p4 > Pa/3 = 2_^ F ai Fp* . (10) 

i 

The transition probability is then the squared modulus of the amplitude (jHJ): 

P a p{L) = 2 [\a aP \ 2 - ReO^a^)] (1 - cos A) + |/3 a/3 | 2 + 2Im(/3^a a/3 ) sin A . (11) 

The probabilities for anti-neutrinos can be obtained from this formula by complex- 
conjugating the parameters a Q/ g and ft a p, i.e. by flipping the sign of the last term in eq. (ITTj) . 
Thus, in the presence of NSI the transition probabilities exhibit CP violation even in the 
one mass scale dominance limit, unlike in the case of pure oscillations. We will exploit this 
property of the considered model in order to reconcile the neutrino and anti-neutrino data in 
short-baseline appearance experiments. Since CP violation comes through the interference 
of a a p and j3 a p terms it is important that both terms be present. It follows from eq. f fTTj) 
that in order to suppress (enhance) the transition probability for neutrinos (anti-neutrinos) 
we need the phase of fi* a a0. a p to be close to 37r/2. Indeed, in our fits to be discussed in the 
following we will obtain numbers close to this value. Note that in general P a p summed over 
either of the two indices will not give unity. This is because NSI break the unitarity of the 
neutrino evolution 48 



B. Application to short-baseline data 

Let us now identify the effective parameters which are relevant for describing the short- 
baseline data used in our analysis, as summarised in tab. HI where we divide the experiments 
into appearance and disappearance searches. 



1. Appearance experiments (/i — > e). In this case we have 

[i& r e4 ' 



o> = F*F% , (12) 





Disappearance 






Appear 


ance 




Experiment 


Ref. 


Channel 


Data 


Experiment 


Ref. 


Channel 


Data 


Bugey 


[4] 


v e -> v e 


60 


LSND 


[2] 


^ ->■ i'e 


11 


Chooz 


M 


v e -)• v e 


1 


KARMEN 


[49] 


^ -> ^e 


9 


Palo Verde 


[50] 


Ve -> ^e 


1 


NOMAD 


f51] 


^ -> ^e 


1 


CDHS 


[6] 


^ ->■ ^ 


15 


MiniB (u) 


[3] 


^ -> ^e 


8 


atmospheric 


[18] 


f M ->■ f/i 


1 


MiniB (i/) 


[1] 


^ "^ p e 


8 



TABLE I: Experiments used in the numerical analysis. The oscillation channel and the number of data 
points used in the fit ("Data") are given. The total number of data points is 115. 
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In deriving this expression from eqs. (fTUj) and (J7J) we have used unitarity of the leptonic 
mixing matrix. Notice that the transition probability P^ e (L) does not vanish in the limit 
L = 0: in that case P„ P = \(3fie\ 2 - This is related to the fact that in the presence of NSI the 



fie 



states \vz) and \vf) are not orthogonal [52], see eq. (J2J). This should lead to a nontrivial 



signal in near detectors in appearance experiments. 



2. Disappearance experiments. We will need the survival probabilities P ee and P MAi . For 
the processes of interest (pion decay, nuclear beta decay and inverse beta processes), both 
neutrino production and detection involve transitions between u and d quarks. For the 
reactor experiments relevant for P ee , production and detection are just inverse processes to 
each other. Therefore the relevant NSI are identical. For P MM the production process is pion 
decay, whereas detection proceeds via neutrino-nucleon capture. Since pions couple only to 
the axial current, vector-like NSI would not contribute to neutrino production but still show 
up in detection. 2 Hence, in this case it is possible that NSI at production and detection 
are different. We will discuss this possibility in appendix \A\ For the moment we simplify 
the discussion and assume that also for v^ disappearance NSI are identical in source and 
detection: 



' fia 



' fia 



(u^ disappearance) . 



(15) 



We would like to thank Joachim Kopp for drawing our attention to this point. 



This corresponds to axial- vector like NSI, e^ = e^ — e^ L \ In this case we find from 
eqs. (HUD and © 

The equality of /3 ee and /3 MM to unity is a reflection of the fact that at L = the survival 
probabilities are equal to one. It should be stressed that this is only correct in the case when 
the neutrino production and detection processes are of the same type, i.e., e^g = s^Bi as m 
the case we consider here. Otherwise the survival probabilities at L = would be smaller 
than one. From eqs. (ITI|) and (fT6l) we find for the survival probabilities 

P^ = 1-2^(1 -0^(1 -cos A), P = e,fi. (17) 

There are no CP violation effects for the survival probabilities, again because we take e£g = 
e^g. Note that eq. (TT71) has the same structure as the disappearance probability for (3+1) 



oscillations with the identification ag — > \Ug^\, see e.g. [53]. 

In appendix |X] we will relax the assumption f[T5"j) and allow e^g ^ s^B f° r ^m disappearance 
experiments. As we will discuss there, assuming e^ g = s^g has a small impact on the global 
fit and therefore we will adopt this simplification for our discussion apart from appendix |A] 
and where stated explicitly. 



As mentioned above, for the disappearance experiments neutrino production and detec- 
tion processes involve transitions between u and d quarks. This holds also for the Mini- 
BooNE and NOMAD appearance experiments. The only exceptions are LSND and KAR- 
MEN. While the detection process still involves u and d quarks (inverse beta decay), the 
production process is purely leptonic (muon decay) 3 , see also eqs. (1A1I) and (1A2I) in the 
appendix. Therefore, all experiments except LSND and KARMEN only depend on semi- 
leptonic NSI, and (under the assumption eq. ffT5|) ) we can drop the distinction between 
source and detector: 

F% = F^ = F£ (semi-leptonic NSI only) . (18) 

This means that the same NSI parameters appear for processes involving a given charged 
lepton a, irrespective of whether the process occurs at neutrino production or detection. 
If this relation holds, we obtain from eqs. (IT2"|) and ffTBl) a factorisation property of the 
parameters a: 

|«Me| 2 = ot e a^. (19) 



3 We not not consider here the decay-in-flight data from LSND, where neutrinos are produced in pion decay. 
Muon decay contributes also to neutrino production in atmospheric neutrinos. We do not expect that 
this will have an important impact on our analysis and neglect this effect in the following, treating all 
NSI effects in atmospheric neutrinos as semi-leptonic. 



Therefore, all experiments except LSND and KARMEN depend on the following set of 
parameters: 

a e , a^ |/3 Me |, 8 = Arg (a>/3* e ) , Am^ . (20) 

For LSND and KARMEN purely leptonic NSI will also contribute, while semi-leptonic 
NSI involving the charged muon, e*^, will contribute to the other experiments but not to 
LSND and KARMEN. This leads effectively to a decoupling of these experiments from the 
others, and we can describe LSND and KARMEN by the three independent parameters 

Kf\,\^\,S LK = Arg(<^r ) (21) 

in addition to the common Am^. Note that in this framework MiniBooNE cannot be 
considered as a direct test of LSND, since due to the different production mechanisms 
different parameters enter the transition probabilities for the two experiments. 

Therefore, the global data in the general (3+1) NSI scenario (under the assumption 
eq. (fl5l) ) depend on the 8 independent parameters shown in eqs. (1201) and (T2TT) . We shall 
denote this case by NSF. In appendix [A] we relax the assumption eq. ffTBT) and argue that 
the effect on the global fit is small. In addition to the general NSF case we will discuss in the 
following also a constrained scenario, denoted by NSI C , which can be obtained by assuming 
that the relation ffT8|) holds for the global data, also including LSND and KARMEN. This 
can be achieved by assuming that all NSI parameters involving the charged muon (leptonic 
as well as semi-leptonic) vanish: 

4 = • (22) 

In this case, we find from eq. ([7]) that F^ = F? = U^ and eq. (TT51) holds for all experiments. 
Note that for the processes relevant to the experiments we consider we always have F? = F®, 
and if eq. (|22|) is fulfilled, also eq. (TT5l) holds trivially. Then the global data depends only 
on the 5 parameters from eq. (1201) and the factorisation (j!9p applies in general. 

The factorisation ffl9l) is analogous to the factorisation property that holds in the pure 
(3+1) oscillations case [53j, and 4|a Me | 2 can be identified with sin 2 26 SBL in the limit /3 Me — ¥ 0, 
#sbl being an effective two-flavour mixing angle. Hence, for /3 Me = we recover the (3+1) 
oscillation case. The fact that the parameter a^ e of the appearance amplitude given in 
eq. f|T9l) is a product of two quantities which are constrained to be small from disappearance 
data leads to the well-known tension between the LSND signal and disappearance data in 
(3+1) schemes. Including now NSI (still subject to the approximation (T2"2"|) ) we obtain two 
more independent parameters, namely |/3 Me | and the phase 5. Since 1/3^1 enters only in the 
appearance probability and is unconstrained by disappearance data, one might expect some 
improvement of the tension between these two data sets. We expect significant improvement 
of the fit in the NSF case, due to the decoupling of the LSND evidence from the rest of the 
global data. 

In the following section II III we will present the results of a fit to data in terms of the 
effective parameters given in eqs. f)20p and (1211) . However in both cases, NSF and NSF, we 
will present also a specific realization in terms of the fundamental parameters U a A and e a p. 
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FIG. 1: Constraint from no-evidence data (NEV) compared to the combined allowed regions from LSND 
and MiniBooNE v data (shaded) at 90% and 99% CL for (3+1) oscillations. We show also the individual 
regions from LSND and MiniBooNE v data. 

III. NUMERICAL RESULTS 



In this section we present the results of fits to short-baseline data in the (3+1) NSI 
framework. The data used in the numerical analysis are summarised in tab. HI In total 
115 data points are used. Technical details on our re-analysis of the experiments can be 
found in 18) and references therein. The constraint from atmospheric neutrino data is 



18 



atmospheric neutrinos 

2 



implemented in the following way. As discussed in detail in 

provide a constraint on the parameter d^ defined in that paper, d^ corresponds to \U^ 
in the (3+1) oscillation case. Using the equivalence of the expressions for the survival 
probability in the (3+1) oscillation and (3+1) NSI schemes in eq. (!T7|) . we identify the 
parameters a M with d^ from 18] and generalise the bound from atmospheric neutrinos now 
to hold for the parameter a M in the presence of NSI. A more accurate analysis of atmospheric 
neutrino data with NSI effects in source and detector is beyond the scope of this work. For 
the MiniBooNE anti-neutrino sample we fit the data shown in fig. 1 of (l|, calibrating our 
prediction to the histogram shown in the lower part of that figure. As we will see later, in our 
model it is not possible to explain the low energy excess seen in MiniBooNE neutrino data 
(and to a lesser extent also in anti-neutrino data) below 475 MeV. Therefore, we follow the 
MiniBooNE collaboration and restrict the analysis to the oscillation-sensitive region above 
475 MeV. 

In fig. [1] we show our fit to MiniBooNE anti- neutrino data compared to the LSND region 
in the (3+1) oscillation framework without NSI. Our result is in good agreement with the 
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FIG. 2: Predicted event spectra at the best fit point to appearance data. Shown are the predictions 
for MiniBooNE neutrino (left), MiniBooNE anti-neutrino (middle), and LSND (right) compared to data. 
Dashed histograms refer to (3+1) oscillations, solid histograms to (3+1) NSI. The dashed vertical lines 
indicate the 475 MeV threshold used for MiniBooNE data. For LSND we show the spectrum with free 
normalisation and the value for the total transition probability (inset). The parameter and \ 2 values are 
given in tab. ILT1 



region obtained in [I]. 4 We observe the overlap with the region indicated by LSND, which 
motivates the combined analysis of the two experiments shown as shaded regions. We also 
show the constraint coming from all the other experiments ("no-evidence" data, NEV), 
which excludes the region to the right of the blue curves. The figure shows the well known 
tension in the (3+1) oscillation fit: the regions touch each other at Ax 2 = 12.7, which 
corresponds to 99.8% CL for 2 dof. 



A. The constrained (3+1) NSI model 

Now we proceed to the (3+1) NSI model, starting with the constrained version NSI C . 
As discussed in section III Bl in this case the global data depends on the 5 parameters from 
eq. (120]) : a e , a^, \f3^ e \, 6, Am^. For appearance data alone the parameters a e and a M enter 
only in the combination \a^ e \ = y/a^a^, while disappearance data do not depend on \Pn e \ 
and S. The limit of (3+1) oscillations is obtained for |/3 Me | = 0. 

First we consider appearance data only, i.e., MiniBooNE neutrino and anti- neutrino data, 
LSND, KARMEN, and NOMAD. The best fit parameters and \ 2 values are given in tab. [Til 



Differences can be attributed to the more sophisticated statistical analysis applied in [1], in contrast to 
the simple A% 2 method based on the Gaussian approximation for 2 dof used here. In particular, Ref. [1| 
obtains a closed contour at 99% CL, while we find only a limit at 99% CL. 
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Data set \ a ^e\ \Pne\ $ Am| x x 2 /dof 

Appearance 0.2075 0.0091 1.5vr 0.1 eV 2 33.5/(37 - 4) 
Global 0.019 0.017 1.3vr 0.89 eV 2 107/(115 - 5) 

TABLE II: (3+1) NSI C best fit parameter and x 2 values for appearance data (LSND, MiniBooNE v, 
MiniBooNE v, KARMEN, NOMAD), and the global data. Corresponding event spectra are shown in figs. 
and[3J respectively. For the global data we have a e — 0.014 and a^ — 0.026 with |a Me | 2 = a e a M . 

and fig.[2]shows the best fit spectra for MiniBooNE and LSND, where the dashed histograms 
correspond to (3+1) oscillations and solid histograms include NSI in addition. The figure 
illustrates the effect of CP violation: the best fit occurs at S = 1.57T (maximal CP violation), 
and we observe the suppression of events for MiniBooNE neutrinos while maintaining the 
signal for anti-neutrinos in MiniBooNE and LSND. Let us mention that KARMEN and 
NOMAD provide an important constraint on the model, in particular they constrain 1/5^1 
to be small. The |/3 Me | 2 term in P Me corresponds to a pure NSI effect, independent of energy 
and distance, see eq. (fTTj) . Therefore the different value of L/E u of these experiments cannot 
be used to circumvent the bounds. This is one reason why a pure NSI explanation is very 
difficult, while in our model we can use the interference terms proportional to Re(a Me /5* e ) 
and Im(a Mej 5* ) in order to circumvent the constraint on the constant term. The relatively 
large value of \a^ e \ ~ 0.2 obtained in the fit is a consequence of the rather small Am^, 
similar to the large mixing angle obtained for an oscillation fit of MiniBooNE anti-neutrino 
data only, see |l|. Such large amplitudes are excluded by disappearance data due to eq. (fl~9l . 
as we will see in the following. 

The last row in tab. [Ill gives the results of the global fit including all data, and the cor- 
responding MiniBooNE and LSND spectra are shown in fig. [3J As mentioned above, disap- 
pearance data constrain |a Me | to be small, and we now obtain |cv^ e | and |/3 Me | of similar order, 
which maximises the CP effect. The phase 5 is now slightly non- maximal, but fig. [3] still 
shows a relevant suppression of the events for MiniBooNE neutrino data, which contributes 
to the overall improvement of the fit. The constraint from disappearance data pushes the 
LSND transition probability to low values, and our best fit prediction Plsnd = 0.19% is 
about 1.8a away from the measured value -Plsnd = (0.264 ± 0.04)% (see inset of right panel 
in fig. [3J. 

In fig. H]we show the allowed regions projected onto the plane of |a Me | and |/3 Me | (left) as 
well as the A% 2 as a function of |/3u e | (right). The improvement of the fit compared to pure 
(3+1) oscillations (1/3^1 = 0) is 

Xmin,(3+l)osc — Xmin,(3+1)NSF = 6-9 (2 dof) , (23) 

where the number of dof corresponds to the two additional parameters of the (3+1) NSI C 
model compared to (3+1) oscillations. Hence, the NSI case is favoured at 97% CL (slightly 
more than 2a) compared to the pure oscillation case. From the right panel of fig. H] we find 
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FIG. 4: Global fit in the (3+1) NSI C framework. Left: allowed regions projected onto the plane of |a Me | 
and |/3 Me | at 90% and 99% CL (2 dof). Right: A% 2 as a function of \/3,j, e \. We minimise over all undisplayed 
parameters. 

that the allowed interval for \f3^ e \ (1 dof) does not include zero at the 2.6a. We conclude 
that (3+1) NSI provides a significantly better fit than (3+1) oscillations. However, despite 
this relative improvement, we stress that some tension remains in the fit also for (3+1) NSI C , 
as discussed below. 
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FIG. 5: Constraint from no-evidence data (NEV) compared to the combined allowed regions from LSND 
and MiniBooNE v data (shaded) at 90% and 99% CL for the (3+1) NSP case. The regions correspond to 
sections of the 4-dimensional allowed regions at fixed values of \fi iie \ and S as given in the legend. These 
values correspond to the ones where the NEV constraint and the LSND+MBP allowed region touch each 
other (at Ax 2 = 11.7). 

Fig. [5] shows the constraint from no-evidence experiments compared to the allowed re- 
gion from LSND combined with MiniBooNE anti-neutrino data. This figure should be 
compared to fig. [T] for (3+1) oscillations. Note that the parameter 4|a Me | 2 used in fig. \5\ 
corresponds to the effective mixing angle sin 2 29$bl used in fig. [1] in the oscillation limit 
|/$ue| — *" 0. The regions shown in fig. Oare sections of the 4 dimensional volume (in the space 
I^Vel, l/^uel) £ 5 Am^) at the fixed values of |/3 Me | and 6 for which the two regions start to touch 
each other. This happens at a A% 2 = 11.7, corresponding to 98% CL for 4 dof. Therefore, 
we observe an overlap of the 99% CL regions. Clearly, while some tension remains in the 
fit, introducing NSI reduces the disagreement between evidence and no-evidence data. 



Let us now address the question of how to realize the effective parameters from eq. (120]) 
in terms of the fundamental parameters U a i and e a p while respecting the phenomenological 
bounds. In fact, it is enough to assume only one single NSI parameter to be different from 
zero, namely e^r. Then, we neglect the quadratic term in the normalisation factor Cg and 
the term containing the product of two small quantities, e^ and U^, in the second equation 
in (TT3|) . As a result, we arrive at the following identification of the parameters relevant to 
our calculation: 



tt M e — U eA U ^ , 



A 



fie 



'efi 



a,. 



\U, 



el 



Oti, 



\U, 



i-i'U 



(24) 
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From tab. [IT] we obtain 

l<l~ l/y ~ 0.017 (25) 

at the global best fit point. This value is in safe agreement with the current bound on this 

Note 



35 



parameter from CKM unitarity and lepton universality, which is of order 0.04 
that the NOMAD experiment sets an even stronger bound of about 0.026 on this parameter. 
Since NOMAD data is explicitly included in the fit it is clear that our result respects even 
this stronger bound. While eq. f)24p is just one simple example, we do not rule out here the 
possibility that other combinations of U^ and e a p (including (3 = s for the sterile neutrino) 
may lead to a similar fit. 

One question that is still to be answered is if it is sufficient to have only one non- vanishing 
parameter e a /3 in order to have physically observable CP violation in neutrino transition 
probabilities. CP violation enters our formulas through the parameter Im (U e /JJ* A e™fy . It 
is easy to see that the phase of the expression in the brackets cannot be rotated away by a 
rephasing of the involved fields, so that CP violation is physical in the case we consider. 5 



B. The general (3+1) NSI model 

Let us now move to the general (3+1) NSI model (NSI 9 ), still assuming e s = e D for v^ 
disappearance experiments, see eq. (TT5l) . As discussed in section iHBl now LSND and KAR- 
MEN data are fitted by the new independent parameters |cCg+ |/?ue+ ^ LK ' while all other 
experiments depend on a different set of parameters a e , a M , \(3^ e \, 5, where the factorisation 
\ a tie\ — y/cte&fi holds. The only parameter in common between LSND/KARMEN and the 
rest is Am^. Therefore, we have now 8 parameters in total. 

In tab. IIHI and fig. [6] we show the result of fits to appearance data only and global data. 
In the first case we obtain an excellent description of MiniBooNE and LSND data. We 
can invoke CP violation to reconcile MiniBooNE neutrino and anti-neutrino data, while 
maintaining an excellent fit to LSND data. 6 However, in the second case, at the global 
best fit point, we observe from fig. that the excess in MiniBooNE anti-neutrino data is 
not explained. The reason is that the constraint coming from MiniBooNE neutrino data 



5 Indeed, U e 4 can be made real, e.g., by rephasing 1/4, and the phase of U^i can be rotated away by rephasing 
\x. However, the parameter e"^ cannot be made real after that because the phases of the e and Vy, fields 
are already fixed. Note that a re-phasing of /i must be accompanied by the corresponding re-phasing of 
Vp due to the standard CC interaction term. One cannot eliminate the phase of e^ff by rephasing the 
fields of u or d quarks that enter into eq. (TTJ) because these fields enter similarly into both standard and 
non-standard interactions, and so the relative phase between these two kinds of terms is unaffected by 
rephasing u or d. 

6 It turns out that allowing for NSI also leads to an improvement in the joint LSND + KARMEN fit. 
Due to the interference terms between c+k anc j ^lk th ere j s more flexibility in the energy dependence of 
the transition probability, which improves the compatibility of the two experiments compared to a pure 
oscillation fit by several units in \ 2 ■ 
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Data set 



K> 



LKi |«LK| rLK 



n 



/xe 



IA 



A'e 



5 



Am| x 



X 2 /dof 



Appearance 0.31 0.029 0.49vr 0.15 0.011 1.5vr 0.13 eV 2 29.4/(37 - 7) 
Global 0.053 0.036 0.39vr 0.010 0.013 1.2tt 0.89 eV 2 95.4/(115 - 8) 

TABLE III: (3+1) NSF best fit parameter and \ 2 values for appearance data (LSND, MiniBooNE v, 
MiniBooNE v, KARMEN, NOMAD), and the global data. Corresponding event spectra are shown in figs.[5J 
For the global data we have a e = 0.010 and a^ = 0.011 with \a UP ] 2 = a P a,, 
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FIG. 6: Predicted event spectra at the (3+1) NSF best fit point to appearance data (dashed) and global 
data (thick-solid). Shown are predictions for MiniBooNE neutrino (left), MiniBooNE anti-neutrino (middle), 
and LSND (right) compared to data. The thin-solid (green) histograms for MiniBooNE show the global fit 
result without the assumption e s = e D for v^ disappearance data (see appendix [XJ. The dashed vertical 
lines indicate the 475 MeV threshold used for MiniBooNE data. For LSND we show the spectrum with free 
normalisation and the value for the total transition probability (inset). The corresponding parameter values 
are given in tab. IIIII 

as well as from the disappearance experiments is much stronger than the positive signal 
in MiniBooNE anti-neutrino data. This is different in the NSF case (c.f., fig. |3]), where 
the signal for MiniBooNE anti-neutrino data is directly linked to the LSND signal which is 
statistically much stronger. Let us mention that should the MiniBooNE anti-neutrino signal 
become more significant in the future, a mechanism similar to the one in the (3+1) NSF 
case (CP violation due to NSI) can be invoked to explain the excess while satisfying the 
bounds. 

In the global fit we find a minimum \ 2 value of 95.4, which corresponds to 



2 _ 2 

Xmin,(3+l)osc Xmin,(3+1)NSI 9 



18.5 



(5 dof) 



(26) 



where the number of dof corresponds to the additional 5 new parameters when extending 
the (3+1) oscillation scheme to NSP. The A% 2 value corresponds to 99.76% CL. Hence, 
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FIG. 7: Global fit in the (3+1) NSI 9 framework. Left: allowed regions projected onto the plane of a^ and 
0\^ relevant for LSND and KARMEN, right: allowed regions projected onto the plane of \a^ e \ and \j3^, e \ 
relevant for all the other experiments. Regions are shown at 90% and 99% CL (2 dof). We minimise over 
all undisplayed parameters. The stars indicate the global best fit point, whereas the triangles correspond to 
the example from eq. 



(3+1) oscillation can be excluded at the 3<r level compared to the NSI 9 case. In contrast 
to the NSI C model, here the tension between appearance and disappearance experiments is 
significantly relaxed, since LSND (which provides the main appearance signal) is decoupled 
from the disappearance experiments. We will quantify this in the next section. Let us 
also mention that there is a nearly degenerate minimum to the one given in tab. IIHI at 
Am| x « 1.8 eV 2 , see fig. [9] below. 

We see from tab. IIHI that the global NSI 9 best fit point requires rather large values for 
the parameters \ct^\ and \(3}^\ relevant for LSND and KARMEN. However, the allowed 
regions shown in fig. [7] extend to rather small values even at the 90% CL. In the following 
we provide one possible realization in terms of fundamental mixing and NSI parameters. It 
turns out that the main difficulty in finding e parameters within the present bounds is to 
obtain \o^\ ^ 0.02 while maintaining \a^ e \ < 0.02, compare with fig. [7J 

Let us take the following e to be non-zero: 



' fis I 



0.05, 



'e/i\ 



0.011 



' flSI 



0.03, 



' fie \ 



0.01 



(27) 



Here the superscript ev indicates the purely leptonic NSI relevant for the muon decay. These 
values are in agreement with the bounds derived in [35|. Using the fit results and eqs. (fT2"|) . 
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dH, dUD we obtain 

Ke\ « (IM -K d s\)\Ue4\™ 0.018, 

|C/ e4 | » 0.116, |£/ M | ^0.205, 

l/VI « |<|~ 0.011, 

K K | « (1*7*1 + l^l)l^4| ^0.027, 

l#fl « Kl + l<el « 0.021 . (28) 

Here we have assumed \U 8 4\ ~ 1 and neglected terms that are quadratic in small quantities. 
Note that we make use of the freedom to choose the phases of the e to suppress |a Me | and 
enhance \afcf\- The x 2 & t this point (with Am^ = 0.98 eV 2 ) is x 2 — 101.0, about 5.6 
units larger than the best fit point from tab. 1111} obtained in the phenomenological analysis 
without taking into account constraints on NSI parameters. Considering that the model has 
8 parameters, a A% 2 = 5.6 corresponds to 69% CL. Therefore, the point from eq. (|28~|) is 
located close to the la volume in the 8-dimensional parameter space. 7 

Let us stress that this point is not optimised and should just serve as an example. Ideally 
a full analysis should be performed in terms of the fundamental parameters U a ± and e a p, 
taking into account the constraints on the latter. We leave such a fit for future work and 
proceed here with the phenomenological analysis in terms of the a and (3 parameters. 

Let us comment now on the relevance of the assumption eq. ffT5|) . of equal NSI at produc- 
tion and detection for v^ disappearance experiments, which we so far have always adopted. 
This assumption can be relaxed by using the fact that vector-like NSI do not contribute to 
pion decay. As discussed in appendix |XJ this allows one to decouple the v^ disappearance 
data. The results of such a fit are shown as the thin-solid (green) histogram in fig. [6J We find 
that the predicted spectra for MiniBooNE are qualitatively very similar to our default NSI 9 
fit. We obtain at the best fit point x 2 — 92.7, to be compared with 95.4 for the standard 
NSI 9 including the assumption eq. (115)) . Hence, the improvement of the fit by relaxing this 
assumption is not significant. The reason is that even with the assumption e s = e D for v^ 
disappearance experiments, the constraints from CDHS and atmospheric data are satisfied 
in the global fit. Therefore, relaxing this assumption leads only to an insignificant improve- 
ment of the fit. Furthermore, as discussed in appendix |XJ this solution requires relatively 
large NSI parameters (of order 0.1) and cancellations between e's and elements of the mixing 
matrix. For these reasons we do not consider this possibility further and stick to the NSI 9 
scenario with the assumption eq. (fl5"|) . 



7 The value of \e~ | given in eq. (f!?T)) is somewhat large (though in agreement with the bounds from 35l|). 
It is used to partially cancel the value of |£/ M 4| in |a Me |. A solution with |e™| = leads to a fit of similar 
quality with Ax 2 ~ 7.7. 
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\U ei U^\ Am 2 41 \U e5 U^\ Am 2 51 6 x 2 /dof 

Appearance 0.397 0.94 eV 2 0.375 1.0 eV 2 l.Olvr 26.3/(37 - 5) 

Global data 0.10 0.15 0.47 eV 2 0.13 0.17 0.89 eV 2 1.69tt 109/(115 - 7) 

TABLE IV: Parameter and \ 2 values of the best fit points in the (3+2) oscillation scheme for appearance 
data from LSND, MiniBooNE v and v, KARMEN, NOMAD (upper part), and global data (lower part). 
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FIG. 8: Predicted event spectra at the (3+2) oscillation best fit point to appearance data (dashed) 
and global data (solid). Shown are predictions for MiniBooNE neutrino (left), MiniBooNE anti-neutrino 
(middle), and LSND (right) compared to data. The dashed vertical lines indicate the 475 MeV threshold 
used for MiniBooNE data. For LSND we show the spectrum with free normalisation and the value for the 
total transition probability (inset). The corresponding parameter values are given in tab. IIVI 



IV. COMPARISON OF (3+1) NSI WITH (3+2) OSCILLATIONS 



In this section we present an update of the (3+2) oscillation analysis from 18|] with 
respect to the anti-neutrino data from MiniBooNE. In this model appearance and global 
data depend on 5 and 7 parameters, respectively, as given in tab. IIVI One physical complex 
phase allows for CP violation in v^ — > v e oscillations 17] and is given by 



5 = arg (u: 4 U, 4 U e ,U; 5 ) 



(29) 



The best fit points for appearance and global data are given in tab. IIVI and the correspond- 
ing event spectra for MiniBooNE and LSND are shown in fig. [HI Considering appearance 
data only, a very good fit is obtained for neutrino and anti-neutrino data thanks to CP 
violation and even the MiniBooNE v and v low energy data can be explained due to some 
cancellations between terms involving the two mass-squared differences and 5 ~ it, as dis- 



cussed in 



18 



However, rather large values of the amplitudes [UeJJ^ and [Ue^U^ are 



18 




FIG. 9: x 2 of global data as a function of Am 2 41 for the (3+1) oscillation, (3+1) NSI C , (3+1) NSP, 
and (3+2) oscillation models. In each case we minimise with respect to all parameters except Am^. The 
horizontal dashed lines indicate the global minima. 



needed, in disagreement with constraints from disappearance data. 

In the global analysis some tension is visible in the spectra shown in fig. |HJ The predicted 
event rate for MiniBooNE neutrino data is somewhat too high and the probability for LSND 
is somewhat low: Plsnd = 0.18%, similar to the (3+1) NSI C case. Furthermore, it is 
not possible to explain the low energy data in MiniBooNE below 475 MeV. The relative 
improvement of the fit due to introducing the second sterile neutrino is 



2 _ 2 

A.min,(3+l)osc A.min, (3+2)osc 



5.0 (4 dof) , 



(30) 



where the number of dof corresponds to the additional 4 parameters introduced by moving 
from (3+1) to (3+2). The value given in eq. (1301) corresponds to 71% CL, to be compared 
to 97% CL for NSP (eq. (JUD) and 99.76% CL for NSP (eq. fl2SD). We conclude that both 
versions of the NSI model (but especially the unconstrained case) offer significantly more 
improvement of the global fit compared to the (3+2) oscillation scheme, when considering 
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without Mini- 



the x 2 g arn P er new parameter. The previous (3+2) value obtained in 
BooNE anti-neutrino data was 81% CL (Ax 2 = 6.1). Therefore, we find that the relative 
improvement of the fit when moving from (3+1) to (3+2) oscillations is smaller with the 
new global data than without MiniBooNE v data. 

The x 2 profiles as a function of Am^ for the four models — (3+1) oscillation, (3+1) NSP, 
(3+1) NSP, and (3+2) oscillation — are shown in fig. [9j In this plot we minimise with 
respect to all parameters except Am^. Note that the maximal x 2 for (3+2) is given by the 
X 2 minimum in (3+1), since the (3+1) solution with Am^ is always available for any value 
of Am^. Note also that there are two degenerate minima for (3+2) at Am^ = 0.47 eV 2 
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(3+1) NSP (3+1) NSP (3+2) oscillations 



Evid. vs no-evid. 
Appear, vs disapp. 



X PG / dof 


PG prob. 


x PG / dof 


PG prob. 


x PG / dof 


PG prob. 


23.3/4 


1.1 x 10~ 4 






26.9/5 


6 x 10~ 5 


11.5/2 


3 x 10~ 3 


3.8/2 


15% 


21.7/4 


2.3 x 10~ 4 


Ax7dof 


CL 


A X 2 /dof 


CL 


AxVdof 


CL 


6.9/2 


97% 


18.5/5 


99.76% 


5.0/4 


71% 



Fit wrt (3+1) osc. 



TABLE V: Comparison of (3+1) NSI and (3+2) oscillations. We show the compatibility of LSND + 
MiniBooNE anti-neutrino data ("Evid.") compared to the rest of the global data ("no-evid."), and the 
appearance versus disappearance experiments. Xpg> number of dof, and the corresponding probability are 
given. The lower part of the table shows the improvement with respect to the (3+1) pure oscillation case. 
We give the improvement in % 2 , where the dof corresponds to the number of additional parameters. 

and 0.89 eV 2 corresponding just to a relabeling of the mass-squared differences, compare 
with tab. IIVI The relative improvement of the three models NSI C , NSI 5 , and (3+2) with 
respect to the (3+1) oscillation case is summarised in the lower part of tab. IVl 

Further insight in the quality of the fit can be obtained by evaluatin g th e compatibility 
of different data sets with the so-called parameter goodness-of-fit (PG) [54]. It is based on 
the x 2 function x PG = xL,min - J2i X?,min 5 wnere xL.min is the X 2 minimum of all data sets 
combined and x\ mm * s ^ ne minimum of the data set i. This x 2 function measures the "price" 
one has to pay by the combination of the data sets compared to fitting them independently. 
It should be evaluated for the number of dof corresponding to the number of parameters in 



common to the data sets, see [54J for a precise definition. 

In tab. ED we show the results of such an analysis, testing the compatibility of evidence 
versus no-evidence data (similar as shown in fig. [5J and appearance versus disappearance 
data. These results indicate that in (3+2) as well as in (3+1) NSI C significant tension 
remains between various data sets. For the general NSI model NSI 9 , however, we find 
excellent agreement between appearance and disappearance data. In this case LSND is 
decoupled from the disappearance experiments and this tension is completely resolved. The 
evidence in MiniBooNE anti-neutrino data is not strong enough yet to show up as significant 
tension in the PG test. We do not perform the evidence versus no-evidence test for NSP, 
since it makes no sense to add LSND and MiniBooNE anti-neutrino data here because, as 
we have seen in fig. EJ the MiniBooNE v excess is not explained in the global fit in this 
model. 

V. SUMMARY AND DISCUSSION 

Recent MiniBooNE anti-neutrino data indicate an excess of v e events, in agreement with 
the LSND evidence for v^ — \ v e transitions. It is known that oscillations with a single sterile 
neutrino at the eV scale are not sufficient to explain the global data. We have investigated 
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the possibility that in addition to a sterile neutrino there are some non-standard neutrino 
interactions (NSI), beyond the Standard Model weak interactions. Since matter effects 
are tiny for the short baselines relevant here, we considered charged-current type NSI in 
the neutrino source and detector, parametrised by e a/ 3 as defined in eq. ([T]). Thanks to the 
interference between NSI effects and oscillations with Am| x ~ 1 eV 2 we obtain CP violation, 
even in the presence of only one mass scale. This effect is used to reconcile the indication 
for v^ — > v e in anti-neutrino experiments (LSND and MiniBooNE) with the absence of a 
signal in MiniBooNE neutrino data. 

We have presented a general parameterisation of the relevant transition and survival 
probabilities in the presence of oscillations (within the one-mass scale approximation) and 
NSI, and we have identified particular combinations of mixing matrix elements U a 4 and 
NSI parameters e a p entering in the probabilities. This drastically reduces the number of 
independent parameters and allows us to perform a general fit to global short-baseline data. 

We have considered two versions of the (3+1) NSI model. In the general case (denoted 
NSI 9 ) we make use of the fact that the neutrino production mechanism in LSND (and 
in KARMEN) is muon decay (purely leptonic), whereas in all other experiments neutrino 
production and detection are semi-leptonic, involving transitions between u and d quarks. 
Therefore, in the presence of suitable NSI parameters we can decouple the transition proba- 
bilities in LSND and KARMEN from the rest of the data. In this case we obtain an excellent 
fit to the global data and the tension between appearance and disappearance experiments 
is resolved. Let us mention that in this case MiniBooNE does not provide a direct test of 
LSND, since different combinations of parameters are relevant for them. Also, in the global 
fit the excess observed in MiniBooNE anti-neutrino data is not reproduced. 

For the second version of the (3+1) NSI model we adopt the assumption that NSI involv- 
ing the charged muon can be neglected. In this case exactly the same NSI parameters are 
relevant for LSND and KARMEN as for all other experiments. In this constrained model 
(NSI C ) we make use of the CP violation due to NSI-oscillation interference to reconcile neu- 
trino and anti-neutrino data. We have shown that in the NSI C model there is a factorisation 
between appearance and disappearance amplitudes, similar to that in the (3+1) oscillation 
scheme. Therefore, it is more difficult to satisfy constraints from disappearance experiments 
and some tension is left in the fit. However, also this model provides significant improvement 
of the global fit compared to the pure oscillation case. 

We have presented the results of our fits in terms of effective parameters, representing the 
specific combinations of NSI parameters entering in the transition probabilities. However, 
for both cases, NSI C and NSF, we have provided also examples of how to realise the required 
parameters in terms of the fundamental mixing and NSI parameters. We have shown that 
values in safe agreement with bounds on the various e's can be found to realise our fits. The 
examples given in eqs. (125]) and f f27|) require e's of order a few xlO -2 . 

We have compared the quality of the (3+1) NSI fits to an updated fit in the (3+2) 
oscillation scheme, which also allows for CP violation due to the presence of two relevant 
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mass scales. Similarly to (3+1) NSI, in (3+2) the appearance experiments can be described 
very well. However, we confirm previous results that for (3+2) oscillations significant ten- 
sion remains in the global fit between appearance and disappearance experiments. The 
improvement of (3+2) compared to (3+1) is not significant, in terms of x 2 gain per new 
parameter. Let us mention also that in none of the scenarios considered here we can ex- 
plain the MiniBooNE low energy excess of events when disappearance data are taken into 
account. Therefore, we follow the strategy of the MiniBooNE collaboration and exclude the 
data below 475 MeV from the analysis, relying on a separate explanation for this anomaly. 
The predictions of our model for future experiments depend on the detailed realization 
in terms of mixing and NSI parameters. In general one may expect some signals in searches 
for deviations from the standard three-flavour oscillation picture in both respects, sterile 
neutrino oscillations as well as NSI. Several proposals to search for sterile neutrinos at the 



eV scale have been presented recently, see for example 55H60J . In 6JJ implications of sterile 



neutrinos for latest cosmological data have been investigated. Recent studies on NSI in the 



context of upcoming and far future experiments can be found, e.g., in [62|-|65|. 
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A sp ecific prediction of our scenario are zero-distance effects in appearance searches 48, 
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since our solutions all include some non- vanishing value of the parameter |^ e |- 
This parameter induces a non-zero transition probability even at zero distance from the 
neutrino source, see eq. ( TTTj) . Hence, the observation of an energy independent appearance 
probability at very short distances is a characteristic signature from this kind of models. 



The idea presented in [60j could be particularly useful to search for this effect, since it would 
allow to map out the E v and L dependence of a v e appearance signal. 

Our model may also provide a signature at the LHC Typically, realising CC-like interac- 
tions as the ones from eq. ([1]) require a charged particle as mediator. The NSI parameters 
e measure the strength of the new interactions relative to the standard weak interaction 
strength set by Gp- Therefore, from our fit results, e ~ 0.01, one expects that the mass 
of a mediator for a dimension-6 operator should be roughly one order of magnitude larger 
than the W boson mass. Hence, one might expect charged particles to show up at the TeV 
scale, with good prospects to be observed at LHC. Let us mention, however, that the results 



of [36j, l37[ suggest that NSI at the level of 0.01 are difficult to obtain from dimension-6 
operators without being in conflict with bounds on charged-lepton processes. As discussed 
there, a possibility to obtain such large NSI would be to go to dimension-8 operators and 
allow for some fine tuning. 
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Appendix A: e 7^ e in i/„ disappearance experiments 

As discussed in section III Bl the fact that pions couple only to the axial-vector current 
may lead to the situation that different NSI contribute at neutrino production and detection 
in i/p disappearance experiments. In particular, vector-like NSI, e^ = £^ R ' + e^ L \ will only 
contribute in the detection process but not at production. This effectively introduces new 
independent parameters which decouple also the v^ disappearance experiments. 

Let us introduce the following abbreviations for the various production and detection 
processes: 

fj, : muon decay 

7r : pion decay . . 

Ne : neutrino-nucleus CC interaction involving an electron 
Nfj, : neutrino-nucleus CC interaction involving a muon 

Relaxing now the assumption eq. ( JT5|) we have the following set of parameters relevant for 
the various experiments (in addition to the common Am^): 



LSND/KARMEN: 


(*fj,e ~ 


_ pM zpNe* 
r fi4 T e4 > 


Pfie ' 


_ V^ ipt* rpNe* 
~~ 2—ii fii ei 


MiniBooNE/NOMAD: 


&[ie ~ 


r [i4 r e4 1 


Pre'' 


_ V^ rp-K rpNe* 
2-^ii fii ei 


reactor: 


^ee — 


1 T?Ne\2 
~ \ r e4 1 ' 


Pee ' 


= 1 


CDHS/atmospheric: 


C/x/i " 




PHH 


__ V^ rpir rpNfi* 
Z_/i fii fii 



(A2) 



where the F a i are defined in eq. (J7J). Now it is possible to have P^ = 1 while allowing for 
a non-zero transition probability in MiniBooNE, which however, requires some cancellation 
between NSI parameters and elements of the mixing matrix. For example, one can take 

</ « UU + <sU; s « , e; s « . (A3) 

This implies a MA , ~ and /3 W ~ 1 and therefore P w ~ 1, as required by the data from 
CDHS and atmospheric neutrinos. On the other hand, F\ « U* 4 and we can have P^ e > 
for MiniBooNE, including the possibility of CP violation. 

The results of such a fit for MiniBooNE are shown as thin-solid (green) histograms in 
fig. El which are qualitatively very similar to our default NSI 9 fit. Reactor experiments still 
constrain the value of |-F e ^ e | to be small, which excludes small values of Am^ where a better 
fit to the MiniBooNE spectrum would be possible (such as for example for the appearance 
data only fit shown in fig. [6]), and we find Am^ ~ 0.9 eV 2 at the best fit point. The spectral 
shape of the signal for such values of Am^ does not allow for a better fit of MiniBooNE 
data even without the constraint from v^ disappearance. Decoupling the v^ disappearance 
data by setting P^ = 1, we find a best fit point with \ 2 — 92.7, to be compared with 

23 



95.4 for the standard NSP including the assumption eq. ( TT5|) . Hence, the improvement 
of the fit by relaxing this assumption is not significant. The reason is that, despite this 
assumption, there is already very good agreement between appearance and disappearance 
data, as discussed in section [TV] Even with the assumption e s = e D for v^ disappearance 
experiments, the constraints from CDHS and atmospheric data are satisfied in the global 
fit. Therefore, relaxing this assumption leads only to an insignificant improvement of the 
fit. 

Furthermore, as mentioned above, decoupling v^ disappearance data requires some un- 
pleasant cancellation. At the best fit point shown in fig. |H]we find U^ ~ 0.26. Therefore, 
to cancel this in F ^ one needs e^f of the same order. Note that the constraints on e^ a 



at the level of a few percent come from pion decay processes [35J, which only apply to 
axial-vector-like NSI, while here we need vector-like NSI precisely to avoid the contribution 
to pion decay. Therefore, such large NSI might be phenomenologically viable, though still 
uncomfortably large. Moreover, there is no reason why they should have values such that 
they cancel against U^, since these are two unrelated quantities. Hence, together with the 
observation that the improvement of the fit is not significant, this motivates us to stick in 
the main text to the assumption eq. (fi~5l) in order to simplify the analysis. 

Appendix B: MINOS results and NSI 

Let us comment briefly on the possibility to apply NSI in source and detector as an 
explanation for recent MINOS results, which indicate a slight difference between neutrino 



and anti-neutrino data [471 ]. where we restrict ourselves here to an effective two-flavour 
framework in the /x — r sector. The mentioned results from MINOS are based on data from 
v^ and v^ disappearance searches. Our formalism allows for CP violation in the survival 
probability provided NSI are different at neutrino production and detection. Hence, one 
can use the same mechanism as discussed in appendix |A] Since pions couple only to the 
axial- vector current, vector-like NSI will contribute only at neutrino detection but not at 



the neutrino source, which is pion decay in MINOS [66]. Note that the one mass scale 
dominance approximation used in this work applies approximately also for MINOS, though 
for Am^. We find that relatively large values of e are needed, of order 0.1, while at the 
same time the combined fit of MINOS neutrino and anti-neutrino data does improve only 
by about 2.5 units in \ 2 - 
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We have also considered the possibility to use NC-like NSI in the context of MINOS 
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, which would induce a non-standard matter effect in the fi — r sector, and which therefore 
could lead to a difference between neutrino and anti-neutrino results. We find that (i) 
diagonal NSI e^,e^ <• 0.5 have a negligibly small effect on the fit, and (ii) off-diagonal 
NSI eJJJ? of the order 0.2 lead to an improvement of the combined fit of about 2.5 units 
in x 2 - These values for e®9, eJJP, e^ are about one order of magnitude larger than the 



bounds from atmospheric neutrino data 70]. We conclude that neither CC nor NC type 
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NSI in the \x — r sector provide a viable explanation for a possible deviation of neutrino 
and anti-neutrino results from present MINOS data. The question whether this conclusions 
holds also in a more general three-flavour (or maybe four-flavour) framework is left for future 



work, see also 71, 72 
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